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Abstract 

In this paper we consider and compare special classes of static theories of gradient elas- 
ticity, nonlocal elasticity, gradient micropolar elasticity and nonlocal micropolar elasticity 
with only one gradient coefficient. Equilibrium equations are discussed. The relationship 
between the gradient theory and the nonlocal theory is discussed for elasticity as well as 
for micropolar elasticity. Nonsingular solutions for the elastic fields of screw and edge dis- 
locations are given. Both the elastic deformation (distortion, strain, bend-twist) and the 
force and couple stress tensors do not possess any singularity unlike 'classical' theories. 



1 Introduction 

The aim of this paper is to give an overview of a special class of generalizations of elasticity. In 
classical continuum mechanics, the elastic continuum is viewed as a collection of particles which 
have only three translational degrees of freedom. The particles are taken without structure and 
are idealized as point masses. 

The first attempts to modify the theory of elasticity were done by taking rotational degrees 
of freedom in addition to the translational ones. In such media, which are called Cosserat 
(or micropolar) media, every material particle is considered as a rigid volume element. Each 
particle has six independent degrees of freedom: 3 micro-rotational and 3 translational ones of 
the center. Thus, a micropolar continuum can be described by geometrical points to which a 
micro-rotation vector is attached. Accordingly, a particle is identified by its position vector and 
its micro-rotation vector. 



'Corresponding author. E-mail address: lazar@lmm.jussieu.fr (M. Lazar). 
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Figure 1: Overview of generalizations of elasticity 



For the calculation of stresses, strains and distortions produced by defects (dislocations, 
disclinations, cracks) such theories break down near the defect lines because they are not valid 
at very small distances. For instance, the elastic stress field of a dislocation has a 1/r-singularity 
and the elastic couple stress field has a l/r 2 -singularity. 

Thus, generalized continuum theories which are able to explain the material behaviour on 
the nanoscale are needed. Such theories can build bridges between continuum and atomic 
physics. On the other hand, the singularities in the elastic fields of defects should disappear in 
a theory which is also valid in the defect core region. The following question arises: What kind 
of generalized elasticity theories are able to do that? 

An extension of elasticity theory which goes into another direction than the Cosserat theory 
is the theory of nonlocal elasticity. The concept of nonlocal elasticity was originally proposed 
by Kroner and Datta [1,2], Edelen and Eringcn [3], Kunin [4] and some others. Nonlocal 
theories were introduced to explain the material behaviour on the nanoscale (e.g., in the core of 
defects). Such a theory considers the inner structure of materials and takes into account long- 
range (nonlocal) interactions. In a continuum theory, nonlocality arises due to the finite range 
of interaction between material points. The long range of interaction implies that the quantities 
such as stress and energy are functionals depending on the motion of all points of the body. 
The nonlocal kernel, which appears in the nonlocal constitutive relation between the stress and 
strain fields, weights the elastic constants of the material. The kernel is supposed to model 
the details of the atomic interaction. The theory of nonlocal elasticity was successfully applied 
to the calculation of the stresses produced by defects (cracks, dislocations, disclinations). The 
main goal is the elimination of the classical unphysical singularities of the stress fields [5-14]. 
Thus the stress fields are smooth even in the core region. In addition, nonlocal elasticity loads 
to a natural fracture criterion which is based on the maximum stress hypothesis [9,10]. So, 
materials will fracture when the calculated maximum stress reaches the cohesive stress of atomic 
bonds. However, the strain field singularities arc still present in nonlocal elasticity. 

Another extension of the classical theory of elasticity is called strain gradient elasticity and 
is very close to nonlocal elasticity. The physical motivation to introduce gradient theories was 
given by Kroner [1, 15] in the early sixties. The strain gradient theories extend the classical elas- 
ticity with additional strain gradient terms. Due to the gradients, they must contain additional 
material constants with the dimension of a length, and hyperstresses appear. The hyperstress 
tensor is a higher order stress tensor given in terms of strain gradient terms. In particular, 
early theories of isotropic, higher-order gradient, linear elasticity were developed essentially by 
Mindlin [16-18], Green and Rivlin [19, 20]. Such gradient theories contain strain gradient terms 
and no rotation vector. In this way, they are different from theories with rotational degrees of 
freedom (Cosserat theory, micropolar elasticity). Only hyperstresses such as double or triple 
stresses appear in strain gradient theories [17]. Gradient elasticity was used to calculate the 
stress and strain fields produced by dislocations and disclinations [21-23]. The gradient elas- 
ticity solutions have no singularity in both the stress and the strain fields. On the other hand, 
in first gradient elasticity the double stresses still have singularities at the defect line [23]. We 
notice that in the special gradient elasticity used by Gutkin and Aifanis [24-26] the stresses are 
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still singular. They only eliminated the strain singularities. Also, the gradient theories latter 
developed and used by Gutkin and Aifantis [21,22,27] and the gradient theory by Lazar and 
Maugin [23] are slightly different. Gutkin and Aifantis used two different gradient coefficients 
and Lazar and Maugin used only one gradient coefficient. 

Eringen [10, 28] introduced a theory of nonlocal micropolar elasticity (or nonlocal Cosserat 
theory). In this theory the couple stress also has a nonlocal constitutive relation in addition to 
the asymmetric force stress tensor. By means of nonlocal micropolar elasticity, the singularities 
in both the force stresses and couple stresses can be eliminated [29,30]. But, the micropolar 
distortion and wryness still possess singularities. This was one motivation for Lazar and Mau- 
gin to introduce a gradient micropolar elasticity [31,32]. The gradient micropolar elasticity 
solutions eliminate these singularities at the defect line. 

Thus, in this paper we only consider various generalizations of classical linear elasticity and 
not physically nonlinear elasticity. Nevertheless, physically nonlinear elasticity is an alternative 
way to eliminate classical singularities at the dislocation line (see, e.g., [33]). But nonlinear 
elasticity is still scale-free and no internal lengths appear. 

The aim of this paper is to give an overview of the solutions of dislocations in gradient 
elasticity, nonlocal elasticity, gradient micropolar elasticity and nonlocal micropolar elasticity 
of Helmholtz type. Gradient theory of Helmholtz type coincides with a special class of first 
order gradient theory Nonlocal elasticity of Helmholtz type means that the nonlocal kernel is 
the Green function of the two-dimensional Helmholtz equation. We give the results for screw 
and edge dislocations. We consider only gradient theories and nonlocal theories with only one 
gradient coefficient and one nonlocal parameter, respectively. 



2 Generalized Elasticity of Helmholtz type 
2.1 Gradient elasticity of first order 

In first gradient elasticity, the strain energy is assumed to depend only on the elastic strain 
and on the first gradients of it [17, 18, 23, 34] 

W = W(E ij ,d k E ij ), (2.1) 

where in the most general case of linear elasticity - the incompatible elasticity - the elastic 
strain is given by 

Ei^^idjUi + d^-E?. (2.2) 

Here Ui is the displacement vector and E[- denotes the plastic strain. The elastic strain (2.2) 
is the symmetric part of the elastic distortion 

/3y = - /?£, (2.3) 

where (3^ is the plastic distortion. As usual, the gradient of the displacement defines the total 
distortion 

0&=djUi. (2.4) 

With Eq. (2.2) the strain energy (2.1) may be rewritten in terms of gradients of the displacement 
and the plastic strain according to 

W = W(d (l u 3) ,d k d (lUj) ,E[ 3 , d k E?). (2.5) 

In the theory of dislocations (see, e.g., [35]), the dislocation tensor is defined by 

dij = ejkidkPu = -ejkidkfiu (2.6) 
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as incompatibility condition and, on the other hand, the compatibility condition for the total 
distortion reads 

tjkidkffi = 0. (2.7) 

By differentiating Eq. (2.6) we obtain the translational Bianchi identity for a^: 

d jaij = 0. (2.8) 

If the plastic strain and plastic distortion are zero, we have compatible elastic strain and 
distortion which are only given in terms of the gradient of the displacement and, thus, the 
dislocation tensor (2.6) must be zero. Such compatible strain gradient theories were investigated 
by Mindlin [16 18]. Because we are interested in dislocations which cause plastic strain, we 
must use an incompatible strain gradient theory [23]. 

In gradient elasticity of Hclmholtz type, the following expression of the strain energy is 
postulated [23,34] 

W =\ a iJ E ij + \ £2 ( d ^ij) (dkE l3 ) , (2.9) 

where the parameter e is called the gradient coefficient which has the dimension of a length. 
The pertinent stress tensors can be defined by taking the variation of W 

dW 

°« ■= (2 - 10) 

dW 2 

T ijk ■= ~ \ = £ O k u tj (2.11) 
d{d k E ij ) 

with the Hooke law for an isotropic media 

<Jij = XSijEkk + 2/^£'i :/ , (2-12) 

where A and /x are the Lame coefficients. For nonnegative strain energy, W > 0, one has 

3A + 2ii>0, /i>0, e 2 > 0. (2.13) 

The stress tensor cry = Uji is a Cauchy-like stress tensor and Tijh = Tjik is the double stress 
tensor. Consequently, the double stress (2.11) is a simple gradient of the Cauchy stress. Thus, 
it is a higher-order stress tensor. Only s 2 is a non-standard coefficient of the theory. 

Further, for vanishing external body forces, the force equilibrium equation can be obtained 
from the principle of virtual work as (variation with respect to the displacement) [16] 

dj(aij - dkTijk) = (2.14) 

and with (2.11) we find 

(l -e 2 k)d o t3 = 0, (2.15) 

where A denotes the Laplacian. Finally, it is convenient to introduce a quantity, which is called 
the total stress tensor, and is defined as 

<Tij = o-ij - d k T ijk = X5 i: jEkk + Z^Eij - e 2 (XS ij AE k k + 2itA£^-). (2-16) 
With this definition Eq. (2.14) takes the form 

djaij = 0. (2.17) 
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If we rewrite Eq. (2.16), we obtain for every component of the Cauchy stress an inhomogeneous 
Helmholtz equation 



(l-s 2 A)a ij =a ij . (2.18) 

The RHS of (2.18) is given in terms of the total stress tensor. 

Using the inverse of Hooke's law for the stress aij and try, it follows that the elastic strain 
can be determined from the equation 

(l-e 2 A)E ij =E ij , (2.19) 

o 

where Eij is the classical strain tensor. 

Using Eqs. (2.2) and (2.3), we obtain the coupled partial differential equation 

(1 - e 2 A) [d {i u j} - /?£.)] = d {l u 3) Pf ij} , (2.20) 

where Ui denotes the displacement field and /3fj is the plastic distortion in classical defect theory 
(see, e.g., [36]). Thus, if the following equations are satisfied 

(l-e 2 A)/3y =kj, (2.21) 

(l- e 2 A)/?£=fe (2.22) 
the equation for the displacement field, 

(l-e 2 A) Ul = u u (2.23) 

is valid for the incompatible case. Thus, for defects (dislocations, disclinations) the inhomoge- 
neous parts of Eqs. (2.22) and (2.23) are fields with discontinuities. The discontinuity of the 
displacement field of a defect is usually represented by a branch cut in order to make the multi- 
valued part to a single-valued one. Thus this discontinuity is for mathematical convenience 
only. This can be accepted because the displacement and the plastic distortion are not physical 
state variables. Therefore, the displacement and the plastic distortion depend on the choice of 
the branch cut. For that reason, they are not unique. 
From Eqs. (2.6), (2.21) and (2.22) we obtain 

(l-e 2 A)a tf =Oij, (2.24) 

where ctij is the classical dislocation density tensor. For a straight dislocation it has the form 

otij = bi®nj5(x)6{y), (2.25) 

where bi and nj denote the Burgers vector and the direction of the dislocation line, respectively. 
Then is the two-dimensional Green function of (2.24) and reads 

a ^^^-l h ^ n o K o(r/e) (2.26) 

where r = \/x 2 +y 2 and K n is the modified Bessel function of the second kind and n = 0, 1, . . . 
denotes the order of this function. 

The limit e — > is the limit from gradient elasticity theory to classical theory of elasticity. 
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2.2 Nonlocal Elasticity 

The basic equations of linear, isotropic, nonlocal elastic solids, for the static case with vanishing 
body force, are [5, 6, 8, 10] 

dM = 0, (2.27) 
a ij (r)= J v a{r-r')af{r')Av{r'), (2.28) 



a. 



IV 



a -""H'kk ■ ( 2 - 29 ) 

Here e is the classical strain tensor, , and ov, are the classical and nonlocal stress tensors, 
respectively. In addition, A and \i are the Lame constants and a(r) is the 'attenuation function' 
called nonlocal kernel. It is important to note that the local material coefficients in Eqs. (2.29) 
are the same which appear in (2.12). 

Eringen and Ari [7] found for the two-dimensional case that an excellent match of phonon 
dispersion was obtained, in the entire Brillouin zone, by 

a(r - r') = Kq (y/(x - x'f + (y - y'Y/e) , e > 0. (2.30) 

e is called the parameter of nonlocality. It is interesting to note that the nonlocal kernel (2.30) 
is the Green function for the two-dimensional Hclmholtz equation 

(l-e 2 A)a(r) =5(x)5(y). (2.31) 

Thus, we call (2.30) nonlocal kernel of Hclmholtz-typc. In this way, we deduce Eringen's so- 
called nonlocal constitutive relation for a linear, isotropic solid with Green's function (2.30) as 
the nonlocal kernel. This kernel (2.30) has its maximum at r = r' and describes the nonlocal 
interaction. The normalization condition for the nonlocal kernel is given by 



X' 



a{r-r')dv{r') = l. (2.32) 



In the classical limit (e — ► 0), it becomes the two-dimensional Dirac delta function 

lim a(r - r') = 5{x - x')5{y - y'). (2.33) 

In this limit, Eq. (2.28) gives the classical expression. We notice that Eringen [8, 10] found the 
two-dimensional kernel (2.30) by giving the best match with the Born-Karman model of the 
atomic lattice dynamics and the atomistic dispersion curves. He used the choice e = 0.39 for 
the length, e = e a, where a is an internal length (e.g. atomic lattice parameter) and eo is a 
constant appropriate to each material. 

Applying the Hclmholtz operator (1 — e 2 A) to (2.28), we find the following inhomogeneous 
Helmholtz equation 



(l-e 2 A)a ij =af, (2.34) 

where is the stress tensor obtained for the same traction boundary-value problem within 
the 'classical' theory. The factor e has the physical dimension of a length and it, therefore, 
defines an internal characteristic length. 

It is obvious that a 'nonlocal' strain tensor does not appear in nonlocal elasticity. Thus, 
the strain and the displacement fields are the classical ones in Eringen's theory of nonlocal 
elasticity. The same is true for the double stresses which occur in gradient elasticity. If we 
identify = Uij, Eqs. (2.18) and (2.34) coincide and the stress Oij also fulfills djOij = in 
gradient elasticity. Then the total stress tensor is equal to the classical stress tensor. Only, the 
stress tensor ov, is the unknown field which must be found. 
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2.3 Screw Dislocation 



We begin with the simplest case of a defect - the straight screw dislocation. First, we need 
the classical quantities because they act as inhomogeneous parts in gradient elasticity as well 
in nonlocal elasticity. The classical stress is given by 

*" = drP=l tr> (2 - 35) 

with the Prandtl stress function 

f = ^lnr. (2.36) 

The modified stress can be expressed in terms of a new stress function F as follows 

a zip = d r F . (2.37) 
For a screw dislocation Eq. (2.18) takes the form 

(l-e 2 A)F = ^ lnr. (2.38) 

The nonsingular solution of the modified stress function for a screw dislocation is given by 

F=^{lnr + K (r/e)}, (2.39) 

which is a superposition of the Prandtl stress function and a gradient part depending on e. 
Thus, it is the two-dimensional Green function of the following partial differential equation of 
fourth order (Hclmholtz Laplace equation) 

(1 - e 2 A) AF = f ib z 5(x)5(y) . (2.40) 

Because of the fact that (2.39) satisfies the inhomogeneous Helmholtz equation in which the 
Prandtl stress function is the inhomogeneous part, one might call (2.39) the Hclmholtz-modificd 
Prandtl stress function. 

Using Eqs. (2.37) and (2.39), the stress field is given by 

^ = ^i{l-^ l(r / e) }. (2.41) 

Unlike the classical result, it does not possess any singularity at r = 0. It has a maximum 
a zip ~ 0.399/j,b z /[2ire] at about r ~ 1.114e. The plot of the stress versus r/e is shown in 
Fig. 2. The stress (2.41) is also the solution of Eq. (2.34) in nonlocal elasticity. In nonlocal 
elasticity this solution was found by Eringen [8] . It is interesting to notice that a solution of a 
screw dislocation in nonlocal elasticity with a nonlocal kernel which is the Green function of a 
diffusion-like equation was earlier given by Eringen [5]. 

By means of the inverse of Hooke's law, the elastic strain is given by 

E zv = ^ 1 -{l- r 1 K 1 (r/e)}, (2.42) 

which is nonsingular at the dislocation line and has its maximum E zip ~ 0.3996 z /[27re] at 
r ~ 1.114 e. The components of stresses and strains in Cartesian coordinates are given in 
Refs. [21-24,37]. 
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Figure 2: The stress <r zip is given in units of /j,b z /[2ire]. The dashed curve represents the classical 
result. 



2.4 Edge Dislocation 

The components of the classical stress of a straight edge dislocation with Burgers vector along 
the x-direction read in cylindrical coordinates 

a rr = ~d r f + ~2 dp V f, a vv = d^ r f, (7 rtp — —d r ^- d^f^j , a zz = v(o rr + a vv ) , (2.43) 

where the Airy stress function is given by 

o uh 

f= X x siny (rlnr) . (2.44) 

2tt(1 - v 7 v ; 

For the modified stress of an edge dislocation, the stresses are given in terms of a stress function 
/ similar to (2.43) 

a rr = ^d r f + ^d* v f, o w =d 2 rr f, a rip = -d r d v f S j , a zz = v(o rr + a w ). (2.45) 

For an edge dislocation (2.18) has the following form 

(l-g 2 A)/ = - ^-g^ sin^(rlnr), (2.46) 

with the Airy function as inhomogeneous part. The nonsingular solution of (2.46) is given 
by [14] 

f = ~2^T^) ™M rlnr + v( 2 - l Kl{r,£) )}- (2 - 47) 

It is a superposition of the Airy stress function and a gradient part. We may call the stress 
function (2.47) a Hclmholtz-modified Airy stress function since it fulfills (Helmholtz bi-Laplace 
equation) 

(1 - £ 2 A)AA/ = d y 6(x)5(y). (2.48) 
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Thus, it is a Green function of this partial differential equation of sixth order. 

By means of Eqs. (2.45) and (2.47), the stress field of an edge dislocation is calculated as 



>rip 



[ib x sin (f r 4e 



'ipip 



&zz = - 



7r(l — v) r 
and the trace of the stress tensor reads 



2,(1-,) r {' -£+«*■/•>}• < 2 -»> 
<*" S "'(l-I KlH£ )}, (2.52, 



^6 x (l + z/) sin<£ r r -i 

fffcfe = ^ \l--K 1 {re)\. (2.53) 

7r(l - f) r I e J 

The stress is plotted over r/e in Fig. 3. Based on the present solution (2.49)-(2.53) we ob- 
serve that the stress fields arc not singular at r = 0. In fact, they are zero at r = 0. 
The radial extremum values are: a rr ~ — 0.259^6a/[2-7r(l — u)e]smip at r ~ 1.494 e, oy v ~ 
0.259/x6 x /[27r(l - v)e] cos ip at r ~ 1.494 £, cr w ~ -0.546^6 a; /[27r(l - z/)e]sin<£ at r ~ 0.546e, 
er zz ~ — 0.399/i6 x ^/[7r(l — i/)e] sintp at r ~ 1.114 e and cr^ — — 0.399/i&a;(l + i/)/[7r(l — v)e] sint/? 
at r ~ 1.114 e. We note that the stresses (2.49)-(2.53) are also solution of Eq. (2.34) in the 
nonlocal elasticity. The solution of an edge dislocation given by Eringen [6] does not correspond 
to the nonlocal kernel (2.30). It is the solution in nonlocal elasticity with the nonlocal kernel 
which is Green's function of a diffusion-like equation. Therefore, it is different. 
The elastic strain of a straight edge dislocation is given by 



E r 



" >]IIP [(1 -2u)-^ + 2K 2 (r/e) + 2*/ K^r/e)), (2.54) 



47r(l — v) r 
b x cos ip r 4e 2 



E rv = ^{l - -L- + 2K 2 (r/e)}, (2.55) 

5 W = -^^y ^{(1 - 2 -) + ^ - ^(r/e) 2(1 - vf-K^r/e)}. (2.56) 



The dilatation reads 

Again the strain fields (2.54)-(2.57) are not singular at r = 0. They are zero at the dislocation 
line. They have the following radial extremum values [y — 0.3): The strain E rr has two 
extremum values of opposite signs (see Fig. 4a), namely, E rr ~ 0.0106a,/ [4-7r(l — v)e] simp at 
r ~ 0.218e and E rr ~ — 0.0546a;/ [47r(l — v)s]sm(p at r ~ 4.130e. In addition, it is equal 
to zero at r ~ 0.677e. On the other hand, E rip ~ 0.2596a,/ [47r(l — ^)e]cos<p at r ~ 1.494e, 
£ w ~ -0.3086^/ [4tt(1 - u)e]sm<p at r - 0.922e (see Fig. 4b) and E kk ~ -0.399^(1 - 
2v)/[2ir(l — v)e] sintp at r ~ 1.114 e. Unlike the classical solution, i? rr is significantly smaller 
than E w in the dislocation core region. The components of the stress and strain tensors in 
Cartesian coordinates can be found in Refs. [14,21-23,25,37]. 
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Figure 3: The components of stress: (a) a vr and (b) a vv are given in units of fib z /[2n(l — v)e\ 
and (c) a zz is given in units of /i6 z ^/[7r(l — v)e] for (p = 3tt/2. The dashed curves represent the 
classical components. 
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3 Generalized Micropolar Elasticity of Helmholtz type 
3.1 Gradient Micropolar Elasticity of First Order 

In first gradient micropolar elasticity (or gradient Cosserat theory), the strain energy is assumed 
to depend only on the micropolar distortion 7^ , the micropolar bend- twist Kij and on the first 
gradients of them [31,32] 

W = W(~/ i j,K i j,dk'yij,dkKij), (3.1) 

where in the most general case - the incompatible micropolar elasticity - the elastic micropolar 
distortion is given by [31,38] 

lij = d j u i + e ijk<fk - ify (3.2) 

Here (fii is the micro-rotation vector and 7^ denotes the plastic micropolar distortion. In 
addition, the elastic micropolar bend- twist (or wryness) reads [31,38] 

= djifi - Kjj , (3.3) 



11 



where nfj denotes the plastic micropolar bend-twist tensor. The total deformation quantities 
are given by 

7^ = djUi + e ijk (p k , (3.4) 

«5 = d m- ( 3 - 5 ) 

For dislocation and disclination it is necessary to use the incompatible gradient theory because 
both defects cause plastic fields. Using Eqs. (3.2) and (3.3), the strain energy (3.1) might be 
written according to 

W = W(d j u i ,d k d j Ui,(pi,d :i ipi,d k d j ipi,'yf j ,dk'rf j ,Kf j ,d k K^). (3.6) 

Thus, the expression (3.6) contains gradients of the displacement and rotation vectors in addi- 
tion to gradients of the plastic fields. 

By means of the plastic fields and the incompatible elastic ones the dislocation and discli- 
nation densities can be calculated as [31,38]. 

a ij = e jkl(9k7il - eUmKmk) = -tjklipkla ~ ^mi) ' O*-?) 

©*j = €jkid k Ku = —ejkidkKtf. (3.8) 

By differentiating Eq. (3.7) and Eq. (3.8) we obtain the translational and rotational Bianchi 

identities for the dislocation and disclination density tensors in the theory of micropolar elas- 
ticity: 

djOtij - e ijk Q jk = 0, (3.9) 

djSij = 0. (3.10) 

The compatibility conditions are 

ejkl (dklli ~ e«m«mfc) = °> ( 3 - n ) 

e 3 kid k nl = 0. (3.12) 

When disclinations are absent, the disclination density tensor 0^ and the plastic micropolar 
bend-twist tensor must be zero and we have Kij = djifi and jfj = /3^. Then Eq. (3.8) is 
the compatibility condition for Kij. The elastic micropolar distortion (3.2) may be rewritten in 
the form 

7i? = Pij + tijkfk- (3.13) 

Then the dislocation density tensor simplifies to 

otij = ejkidkPii = -tjkldkPu- (3.14) 

In gradient micropolar elasticity of Helmholtz type the strain energy, W, is assumed to be 
given as [31, 32] 

W = i (Titfij + i HijKij + ^£ 2 (d k aij) (d k jij) + ^£ 2 (dkmj) (d k Kij), (3.15) 
where k is the only one gradient coefficient. By means of the strain energy (3.1), the asymmetric 
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stresses are defined by 

_ dW_ 

° 13 :_ 



(3.16) 



dW 

■= ^? (3 - 17) 

dW 2 

T i?fc : = 7777; \ = £ Ok&ij, (6. LIS) 

<j{Okli 3 ) 

\jk ■= 7777!— - \ = £ 2 Sfc^ij- (3-19) 

Here (Xy and Hij are the force stress tensor and couple stress tensor, respectively. The higher 
order stresses r^fc and \jk are called double force stress and double couple stress, respectively. 

The constitutive relations for isotropic micropolar elasticity (see, e.g., [10,39,40]) have the 
form 

°y = A %7fefe + (m + r /)7ij + (M - »y)7jii (3.20) 
Mij = a <%«fcfc + ((3 + + ((3- ^Kji. (3.21) 

Here A, /j,, rj, a, (3 and 7 are the 6 material constants of the isotropic Cosserat continuum. For 
nonnegative strain energy, W > 0, we have 

3A + 2/x>0, ^ > 0, 77 >0, 

3a + 2/3>0, (3>0, 7 >0 (3.22) 

and 

s 2 > 0. (3.23) 

Using the six material constants of a Cosserat continuum (micropolar medium) , two character- 
istic lengths I and h can be defined by [39, 40] 

I- = + ^ + , 2 = « + (3 . 24) 

4/xr/ 4?7 

which are Cosserat intrinsic lengths. The above equations (3.20) and (3.21) may be expressed 
in terms of stresses 

Jij = A' SijVkk + (p' + + (p' - v'fau (3.25) 
Kij = a' 5 t:i ii kk + (/?' + + (J3' - 7')%. (3.26) 

Here 

^=h 2ri= h 2l '=h (3 - 27) 

A ' = _ 2 A1 (3A + 2/i)' a ' = ~2/3(3a + 2/?)- (3 ' 28) 

After variation of (3.1) with respect to the displacement and the rotation vector ipi, the 
following equilibrium conditions follow 

a, (cry - a fe T yfe ) = 0, (3.29) 
dj{nij - dkXijk) - eijk((Tjk - diTjki) = 0. (3.30) 
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We may define the total stress and total couple stress tensors 

crij = <Jij - dknjk, (3.31) 
Mjj = Mij _ dk\jk- (3.32) 
So, Eqs. (3.29) and (3.30) take the form 

djOij = 0, (3.33) 

dj^ij - eijkVjk = (3.34) 

or in terms of cr^ and /Uy 

(1 -s 2 A)d j a ij = 0, (3.35) 

(1 - e 2 A) (djUij - eyfc<7jfc) = 0. (3.36) 

From Eqs. (3.31), (3.32), (3.18) and (3.19) the force stress tensor cry and couple stress tensor 
Hij satisfy the inhomogencous Hclmholtz equations 

(l-e 2 A)a ij =a ij , (3.37) 
(l-£ 2 A)/iy =^j, (3.38) 

where the RHS are given in terms of total stress and couple stress tensors. 

Using Eqs. (3.25) and (3.26), the inverse of Hooke's law for the stress and cry, it follows 
that the elastic strain can be determined from the equation 

(l-e 2 A) 7ij =7y, (3.39) 
(1 -e 2 A)n l0 = Kij, (3.40) 

where jij and Kij are fields in 'classical' micropolar elasticity. These are coupled partial differ- 
ential equations. In fact, substituting Eqs. (3.2) and (3.3) in (3.39) and (3.40), we obtain 

(1 - e 2 A) [diUj + e ijk ipk - = diUj + e ijk ^ k - jfp (3.41) 
(l-e 2 A)[d m -4]=dw j -°4, (3.42) 

where <£j denotes the rotation field, 7^ is the plastic micropolar distortion in defect theory (see, 
e.g., [38]). Thus, if and if only the following equations are fulfilled 



(1- £ 2 A) 7 £ = 7£, (3.43) 
{l-e 2 A)4=°^ (3.44) 

the equations for the displacement and rotation fields, 

(l-e 2 A)u l = u l , (3.45) 
(1- £ 2 A)^ = ^, (3.46) 

are valid for the incompatible case. Since nfj and are the plastic fields of disclinations, they 
are zero for dislocations. In this way, Eq. (3.46 follows directly from (3.42) for dislocations. 
Thus, in dislocation theory only the inhomogeneous parts of Eqs. (3.45) and (3.43) are fields 
with discontinuities. 
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In addition, we obtain for defects in first gradient micropolar elasticity from Eqs. (3.7), 
(3.8), (3.39) and (3.40) 

(l-e 2 A)a t3 =a tf) (3.47) 
(l-e 2 A)e ij = e ij , (3.48) 

o 

where 8y is the classical disclination density tensor. For a straight disclination it reads 

Q ij =n i <»n j 5(x)5{y), (3.49) 

where fij denotes the Frank vector. Consequently, the dislocation density and the disclination 
density tensors of a straight defect have in first gradient micropolar elasticity the following form 

an = ^3 h <g> rij K (r/e) (3.50) 

&ij = Sli <8> nj K (r/e). (3.51) 

In the limit e — > 0, we recover micropolar elasticity. Thus, it is the 'classical' limit in 
gradient micropolar elasticity. 

3.2 Nonlocal Micropolar Elasticity 

The basic equations in the isotropic theory of nonlocal micropolar elasticity are given by [10, 28] 

djaij = 0, (3.52) 
djUij ~ £ijk&jk — 0, (3.53) 

<7«(r)= / a(r-r')alf(r')dv(r'), (3.54) 
Jv 

Aty(r) = J a(r-r')$\r')dv(r'). (3.55) 

Using the nonlocal kernel (2.30) and applying the Hclmholtz operator to (3.54) and (3.55), one 
obtains 

(l-e 2 A)a ij = ag 1} , (3.56) 
(l- £ 2 A)^=^f, (3.57) 



where ajp and n\p are the classical local stress and couple stress tensors given by 

= A ME 5 + (M + Vhf + (/* - ^/)7jf \ (3-58) 

,( cl ) - ^,A..^( cl ) _i_ ffli^wM _L Cfl _„A*-( cl ) 



fell (cl) 

Here 7^ and 7^ are the local micropolar distortion and micropolar bend-twist tensors, 
respectively. It is important to note that the local micropolar material constants in Eqs. (3.58) 
and (3.59) are the same which appear in (3.20) and (3.21). 

Again if wc identify a^p = Oij and jifP = ji-, Eqs. (3.37) and (3.38) coincide with (3.56) 
and (3.57), respectively. In addition the stresses fulfill the equilibrium conditions (3.52) and 
(3.53) in gradient micropolar elasticity, too. No double force and double couple stresses appear 
in nonlocal micropolar elasticity. 
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3.3 Screw Dislocation 

The expressions for the force stresses and couple stresses of a screw dislocations in a Cosserat 
continuum or micropolar medium were given by Kessel [41] (see also [39,40,42,43]). In mi- 
cropolar elasticity, the asymmetric force stresses of a straight screw dislocation may be given 
in terms of the stress function 

& zip = d r F_ , <t vz = d r F + , (3.60) 

o 

where the two stress functions F± are given by 

F± = ^{»lnr±r]K (r/h)}. (3.61) 

They fulfill 

(1 - h 2 A) AF ± = b z {(^ T »?)[! - h 2 A] ± V }S(x)5(y) . (3.62) 

In gradient micropolar elasticity, the asymmetric force stress can be expressed in terms of 
new stress functions F_ and F + : 

o zv = 9 r F_ , a vz = d r F+ . (3.63) 

If we substitute (3.60), (3.61) and (3.63) into the Helmholtz equation for the force stress ten- 
sor (3.37), we obtain for the modified stress functions the following inhomogencous Helmholtz 
equation 

(l-e 2 A)F± = ^-{^\nr± V K (r/h)} 1 (3.64) 

where the inhomogeneous part is given by the stress functions (3.61). The nonsingular solution 
of (3.64) is [31] 

, 7 2 

F ± = ^H lnr + *o(r/e)] ± V i K o(r/h) - K {r/e)]}. (3.65) 

It is a more complicated superposition of the micropolar stress function (3.61) and a gradient 
coefficient depending term that in strain gradient elasticity 1 . In addition, they satisfy 

(l -£ 2 A)(l - h 2 A)AF± = 6 z {(^t^)[1 - h 2 A}±f]}S(x)S(y). (3.66) 

Using (3.63) and (3.65), the asymmetric force stress is found as 

r -K x {r/h)- r -K x {r/e)]}, (3.67) 

LK^r/V-'-K^r/e)}}. (3.68) 

They are zero at r = and have extremum values near the dislocation line. The cxtrcmum val- 
ues depend on h and e. For instance, with h — 2e and \i — 3ry we obtain: a zip ~ 0.499 fib z /[2ne] 
at r ~ 1.024e and a vz ~ 0.303 iib z /[2Tie\ at r ~ 1.312e. Thus, o zv > a vz . The stresses are 
plotted versus r/s in Fig. 5. We note that the stresses (3.67) and (3.68) are also solutions of 
Eq. (3.52) in nonlocal micropolar elasticity. In nonlocal micropolar elasticity these solutions 

1 In [31, 32] we used the notation k = 1/e, r = l/h and £ = l/l. In the current notation the pre-factors and 
the limits to micropolar elasticity are easier to handle. 
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Figure 5: The components of the asymmetric stress of a screw dislocation: (a) a zip and (b) a vz 
are given in units of [ib z /\l'ne\ with h = 2e and fi = 3r/. The dashed curves represent the 
micropolar results. 



were calculated by Povstenko [29] (He used the transposed tensors.). In addition, in the limit 
to micropolar elasticity e — > 0, we recover in (3.67) and (3.68) Kessel's result [41]. 

By means of the inverse of the constitutive relation (3.20), the micropolar distortion is 
obtained as 



h l r 

47r r I 



7 vz = ^-{\l--K 1 (r/e) 
4-7T r L L e 



l--K 1 (r/e) 



h 2 



h 2 



h 2 



LK^r/^-^K^r/e)]} 
LK^r/V-^K^r/s)]} 



(3.69) 
(3.70) 



They are zero at r = and have extremum values near the dislocation line. Again, the 
extremum values depend on h and e. For instance, with h = 2e we have: ^ Zip ~ 0.701 6 z /[47re] 
at r ~ 0.941 e and 7 VZ ~ 0.163 & z /[47re] at r ~ 3.395 e. Thus, it yields j zip > -f vz . In addition, 

7 VZ ~ 0.151 b z /[4ire] at r ~ 2.229 e. The micropolar distortions (3.69) and (3.70) are plotted in 
Fig. 6 versus r/e and with h — 2e. It is interesting to note that the elastic micropolar strain 
7( Z(p ) has the same form as the elastic strain E zip calculated in gradient elasticity (see Eq. 2.42). 
In the limit e — > 0, we recover in (3.69) and (3.70) Minagawa's result [42]. 
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Figure 6: The components of the micropolar distortion for a screw dislocation: (a) j zlp and 
(b) 7<^z are given in units of 6 z /[47re] and h = 2e. The dashed curves represent the micropolar 
results. 



The micro-rotation of a straight screw dislocation reads [31] 



(fir = l V z = -d r F + {n = i,n=i). 



(3.71) 



On the other hand, ip r follows directly from Eq. (3.46) with (p r = "f vz . The micropolar bend- 
twist or wryness is given in terms of (3.71) 

1 



r 



So, we find 



Air r 2 



{> 



1 



h 2 -e 2 



k <p<p = ^{ 1 - 7^72 [ hr ^1 - er K * ( r / £ )] } • 
The trace of the micropolar bend-twist tensor reads 



(3.72) 



hrK^r/h) - srK^r/e) + r 2 (K a (r/h) - K (r/e))]}, (3.73) 



(3.74) 



Kkk — . 



An h 2 - e 2 



K (r/h)-K Q (r/s) 



(3.75) 
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The components K rr , k vv and K kk have a maximum at the dislocation line. In fact, one obtains 
at r = 0: 

1 b z 1 h 

K rr 2 K kk ^ ^ ^ _ ^ In g . ^ ttj 

For the choice h — 2e we have at r = 0: 

1 b 



Kkk ~ 0.116^. (3.77) 

Thus, all 1/r 2 - and lnr-singularities, which are present in the results of Cosserat theory, are 
eliminated. In addition, the bend-twist is modified in the region r/s < 6. The shape of the 
bend-twist in gradient micropolar elasticity looks more physical than the bend-twist in classical 
micropolar elasticity (see Fig. 7). In addition, it can be seen that n rr is negative in the far 
field expression and positive near the dislocation line. The point at which it changes sign is for 
h = 2e: r ~ 3.395 e (see Fig. 7a). In the limit e —> 0, we recover in (3.73) and (3.74) Minagawa's 
result [42]. 

Using Eq. (3.21) one obtains from the micropolar bend- twist tensor the couple-stress tensor. 
The non- vanishing components of the micropolar couple stress tensor are 

Mrr = "fe^ ~ [^Kiir/h) - erK^r/e)) + r 2 (K (r / h) - K (r/e))] }, 

(3.78) 

M W - j^-s [pihrK^r/h) erK^r/e)) | (K (r/h) K (r/e))] }, (3.79) 

K {r/h)-K (r/e)\. (3.80) 



b z a 1 

Hz 



Air h 2 - e 

The trace of the couple stress tensor is given by 

6 z (3a + 2/3) 1 



V- kk 



4tt h 2 



K {r/h)-K Q (r/e) . (3.81) 



All components of the couple stress tensor are nonsingular in the gradient micropolar elasticity. 
In fact, one obtains at r = 0: 

a a a b z 1 h . . 

(J-zz = -5- — = -5- — fJ-rr = nf} , » Mfefc - a — —— — In - . 3.82) 

p + a p + a 2p + 3a 2tt 2{h z — e z ) e 

The micropolar bend-twist and the couple stress tensors of a screw dislocation are symmetric 
tensors. The far-fields of (3.73) and (3.74), and (3.78)-(3.80) agree with the result given by 
Kessel [41] and Minagawa [42]. In the limit e — > 0, we recover in (3.78)-(3.80) Kessel's result [41]. 

The couple stresses (3.78)-(3.80) are the solutions in nonlocal elasticity. They fulfill together 
with the stresses (3.67) and (3.68) the moment equilibrium (3.52) in nonlocal micropolar elas- 
ticity. Using the relation K2 (x) — Kq (x) = 2/x K\ (x) , it can be shown that the couple stresses 
agree with the couple stresses in nonlocal micropolar elasticity calculated by Povstenko [29] 
up to a pre- factor 1/2 which is missing in his formulas. This is because he used Nowacki's 
expressions [39] for the couples stresses of a screw dislocation in Cosserat theory which are the 
source of this missing factor 1/2. The results of Kessel [41] and Minagawa [42] have the correct 
pre-factor. The components of the force and couple stress, distortion and bend-twist tensors in 
Cartesian coordinates are given in Ref. [31]. 

3.4 Edge Dislocation 

The expressions of the force stresses and couple stresses of a straight edge dislocations in 
micropolar elasticity (Cosserat theory) are given by [39-41,44,45]. We use the stress function 
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Figure 7: The components of the micropolar bend- twist tensor for a screw dislocation: (a) n rri 
(b) Ktptp and (c) Kkk are given in units of b z /[Aire 2 ] and with h — 2s. The dashed curves represent 
the classical (micropolar) components. 
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formulation given by [46,47] (see also [38,40,43]), which involves two stress functions / and 'J', 
which are connected with the force stresses and couple stresses by 

a rr = \d r f+± d 2 J - \ d 2 J + 1 d v % (3.83) 

o w = d 2 rr °f + - d 2 ^ - \ djl, (3.84) 

a rv = — d%) + 1 + 9 r 2 r *, (3.85) 

a vr = - \ d%) + 1 - J c? r * - 1 5^*, (3.86) 

» tr = d r *, fi Z(p = l -d^ (3.87) 

o o o ° 

and a zz = v{o rr + a vv ). Here / is the Airy stress function known from classical elasticity 

o 

and \I/ is the stress function introduced by Mindlin [46], Schaefer [47] and Carlson [48]. For a 
straight edge dislocation the so-called Airy-Mindlin stress functions reads [39,40] 

o uh 

f = ~ 2 J l X _ v) smy(r-lnr), (3.88) 



| = _(P + l)b x cos Lp 



2ir 
with 



(l- r - Kl {r/l)) (3.89) 



(l-/»A)* = -^±^^. (3.90) 

With the help of Eq. (3.90) and some simple manipulations, we find the Cauchy-Riemann 
relations 

O r (l-l^ = -^± 1 -d v (Af), (3.91) 

^(l-PA)^ = ^- )dr (Af). (3.92) 

By construction this stress function ansatz fulfills the equilibrium conditions (3.52) and (3.53). 
All components of the force stress tensor have 1 / r-singularities and those of the couple stress 
tensor have l/r 2 -singularities. 

In gradient micropolar elasticity we make the following ansatz 

a rr = l - d r f + 1 d 2 w f - - r d 2 ^ + 1 d v V, (3.93) 

^ = 0+;^-^,$, (3.94) 

<r rv = ~ d 2 v f + ^ 8 v f + 9 2 r *, (3.95) 

a vr = -\d%f + 1 d v f -\d r *-^ d 2 ^, (3.96) 

H tT = dr^, /J,z<p = \ (3.97) 

and a zz = v(o rr + er w ). It has the same form as the stress function ansatz in micropolar 
elasticity. Here / and are the modified stress functions which must be calculated. So, we 
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obtain for the modified stress functions the following inhomogeneous Helmholtz equations 

(l-e 2 A)f = -^-^— ) sin^(rlnr), (3.98) 
(l-^ = -^±*^(l-^ l( ,//)). (3.99) 

The nonsingular solutions are (see [31,32] for technical details) 

/ = -^^T sin *,{r In r + ^ (l - \ K^r/e)) }, (3.100) 



2tt(1 -v) r I r \ e 

(P + j)b x coslp r 1 



2tt 



I 1 -p^[j^W ; ) - Jtfi(r/e)]}. (3-101) 



Again, they are superpositions of the classical stress functions (3.88) and (3.89) and the gradient 
parts which depend on the gradient coefficient e. Therefore, the stress function (3.100) has the 
same form as in strain gradient elasticity (see Eq. (2.45)). The stress function (3.101) also 
satisfies the following bi-Helmholtz equation: 

(1 - I 2 A) (1 - e 2 A)* = -(/? + -i)b x d x S(x)6(y) . (3.102) 

Thus, it is the corresponding Green function. In addition, the stress functions ^ and / are 
related by the following Cauchy-Riemann relations 

d r (l - e 2 A) (1 - l 2 A)V = \ 8,(1 - e 2 A)Af , (3.103) 

\ 3,(1 - e 2 A) (1 - I 2 A)* = d r (l - e 2 A)Af. (3.104) 

Using the Airy-Mindlin stress functions (3.100) and (3.101), we find for the asymmetric 
force stress of a straight edge dislocation 

Hb x sin^ r As 2 \ 

+ !^{2 - Wr/0 - tf a ( r , e)] }, (3.105) 

- ^{l + ^ - 2^(r/ £ ) - 2 I K l{ r,e)\ 

2tt{1 - v) r V r z e > 

- ^ - j^W) - **/.)]}, (3-106) 
i^ 2 ?^- £ + «■<•■/•>} 

(3.107) 



fib x cos y> r 4s 



27r(l — r I r 
(/3 + 7)fr* cosiy9 



f 4e 2 i 
{!- — + 2^(r/e)} 



/i6a;^ sin ip 
7r(l — i/J r 



2 , ,., {2- i ^[^ 2 (r/0-if2(r/ e )]} ) (3.108) 
{l-^i(r/e)}, (3.109) 
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and 



fib x (l + v) sin Lp 

Okk — f. r 

7r(l — v) r 



{l-^#i(r/e)}. (3.110) 



The stress fields have no artificial singularities at the dislocation line and the extremum stress 
occurs at a short distance away from the dislocation line (see Figs. 8). The stresses (3.109) and 
(3.110) are not influenced by the Cosserat constants and, thus, they have the same form as the 
corresponding expressions (2.52) and (2.53) in strain gradient elasticity. In Eqs. (3.105)-(3.108) 
the first parts agree with the stresses (2.49)-(2.51) in strain gradient elasticity and the other 
parts are modified stresses due to the gradient micropolar elasticity. Near the dislocation line 
one can see the difference between the stresses calculated in strain gradient elasticity and in 
gradient micropolar elasticity. The extremum values of these stresses are changed. They are 
higher or lower than the extremum values in strain gradient theory (see Fig. 8). In addition, 
the values and the positions of the extremum values depend strongly on the material constants 
of gradient micropolar elasticity. 

Again, the stresses (3.105)-(3.109) are solutions of Eq. (3.52) in nonlocal micropolar elastic- 
ity. This nonlocal solution coincides with the solution given by Povstenko and Matkovskii [30] . 
In the limit e — > 0, we recover in (3.105)-(3.109) the micropolar result given by [39,41]. 

Using (3.26), we obtain the elastic micropolar distortion from the asymmetric force stress 
tensor 

Irr = —. ; ^{(1 - 21/) - ^ + 2K 2 (r/e) + 2v T - K^r/e)} 

47r(l — v) r I r z e > 



47r(l — v) r I r 2 
; + 7> 

47T/Z 

b x sin ip r 4, 



+ iJL ^ ^ - j£,[*W0 - **/«>]}, (3.111) 



(P + j)b x simp 



4tth r 3 I I 2 — e 
b x cos <p r 4e 



{ 2 -jT=^[ K ^/l)-K 2 (r/e)]}, (3.112) 



47r(l — v) r V r 

((3 + j)b x cos Lp 



{l-^- + 2K 2 (r/s)} (3.113) 



{ 2 - [#»(r/Z) - K 2 (r/e) + 1±t (L Kl (r/l) T - K^r/e))] } 



47r/x r 3 I I 2 — e 2 L 2i] \ I 

b x cos ip f 4e 



47r(l —v) r 

((3 + j)b x cos Lp 



and 



{l-^ + 2K 2 (r/e)} (3.114) 
{ 2 - ]2^2 [ K ^/l) K 2 (r/e) + ^ (y ^(r/l) r - K^r/e))]}, 



ub x (l — 2v) simp r r rr , , ,"i ,„ . 

^ = -^rz^-r{ 1 -e^ r /4 (3 - 115) 

It can be seen that the elastic micropolar strain 7^) is influenced by the Cosserat constants /?, 
7 and I and, therefore, it is different from the elastic strain (2.54)-(2.56) of an edge dislocation 
in strain gradient elasticity. Only the first parts of Eqs. (3.111)— (3.114) agree with the elastic 
strains (2.54)-(2.56) in strain gradient elasticity. The other parts of Eqs. (3.111)— (3.114) are 
modified stresses due to the gradient micropolar elasticity. On the other hand, the micropolar 
dilatation (3.115) has the same form as the dilatation (2.57) in strain gradient elasticity. In the 
limit e — > 0, we recover in (3.111)— (3.114) the micropolar result given by [44]. 
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Figure 8: The components of stress: (a) a rr and (b) a vv are given in units of /z6 z /[27r(l — v)e] for 
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dashed curves the stresses in strain gradient elasticity. 
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Using Eqs. (3.97) and (3.101), we find for the micropolar couple stresses 



^zip 
Htpz 



ifi + l)b x cosy 
2tt r 2 

(/3 + 7) b x siny 

2tt r 2 
/?-7 
/3 + 7 

(3 + 1 



1 



/ 2 -e 2 
1 



/ 2 -e 2 



l^zcp- 



IrK^r/l) - erK^r/e) + r 2 (K (r/l) - K (r/e))] }, 

(3.116) 

IrK^r/l) -£r^i(r/£)]}, (3.117) 

(3.118) 
(3.119) 



Due to the special choice of cylindrical coordinates the couple stresses (3.116)-(3.119) have an 
artifical discontinuity at r = 0. At r = the radial part of the couple stresses is finite and is 
multiplied by cos p or sin p>. In Cartesian coordinates the discontinuity disappears (see [32]). 
The couple stresses (3.116)-(3.119) are in agreement with the nonlocal couple stresses calculated 
by Povstenko and Matkovskii [30]. In the limit e — > 0, we recover in (3.116)-(3.119) the 
micropolar result given by [41]. 

With (3.26) we obtain from the micropolar couple stresses the corresponding micropolar 
bend-twist (wryness) 



bz_ cosy r 
2tt r 2 I 



1 



Z 2 - 



b x sin ip 



IrK^r/l) - erK^r/e) + r 2 (K {r/l) - K (r/e))] }, (3.120) 
IrK^r/l)- erK^r/e) }. (3.121) 



The micropolar rotation can be calculated from the micropolar bend-twist. For example, if we 
use 



= rjn z 



,d<p, 



wc obtain 



•£z = - 



b x cos ip 
27 ~ 



i 1 -JI^[ lrK ^ r / 1 ) - zrK^r/e)] }. 



(3.122) 



(3.123) 



For e — > 0, (3.123) coincides with Nowacki's result [40] calculated in the theory of micropolar 
elasticity. The components of the force and couple stress, distortion and bend-twist tensors in 
Cartesian coordinates can be found in Ref. [32]. 



4 Conclusion 

In this paper, we have given an overview over static theories of gradient elasticity and nonlocal 
elasticity of Helmholtz type. We have discussed the theory of defects (dislocations, discli- 
nations) in (gradient) elasticity as well as (gradient) micropolar elasticity. Such theories are 
incompatible gradient theories because the elastic distortion quantities are not simple gradients 
of the displacement and the rotation. We have shown the equivalence between the nonsingu- 
lar stresses of screw and edge dislocation in gradient elasticity and the corresponding nonlocal 
stresses. In addition, we have investigated the relationship between gradient micropolar elastic- 
ity and nonlocal micropolar elasticity of Helmholtz type. For the static case, we have proven the 
equivalence of the force stresses and couple stresses calculated in gradient micropolar elasticity 
and nonlocal micropolar elasticity. We have considered straight screw and edge dislocation. In 
order to fulfill the equilibrium conditions we have used the stress function method. These stress 
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functions are Green's function of partial differential equations of higher order. The calculated 
force and couple stresses do not possess singularities in the dislocation core region. For the 
gradient theories we have calculated the elastic distortion and bend-twist tensors of screw and 
edge dislocations. Unlike the nonlocal theories, where these fields still possess singularities at 
the dislocation line, the quantities calculated in the gradient theories are nonsingular. All fields 
calculated in the theories of gradient elasticity or gradient micropolar elasticity have the correct 
limits to classical elasticity or to micropolar elasticity. 

Acknowledgement 

This work was supported by the European Network RTN "DEFINO" with contract number 
HPRN-CT-2002-00198. 

References 

[1] Kroner, E., Datta, B.K., 1966. Nichtlokale Elastostatik: Abteitung aus der Gittertheorie. 
Z. Physik 196, 203-211. 

[2] Kroner, E, 1967. Elasticity theory of materials with long range cohesive forces, Int. J. 
Solids Struct. 3, 731-742. 

[3] Eringen, A.C., Edelcn, D.G.B., 1972. On nonlocal elasticity, Int. J. Engng. Sci. 10, 233- 
248. 

[4] Kunin, I. A., 1983. Elastic Media with Micro structure II: Three- Dimensional Models, 
Springer, Berlin. 

[5] Eringen, A.C., 1977. Screw dislocation in nonlocal elasticity, J. Phys. D: Appl. Phys. 10, 
671-678. 

[6] Eringen, A.C., 1977. Edge dislocation in nonlocal elasticity, Int. J. Engng. Sci. 15, 177-183. 

[7] Ari, N., Eringen, A.C., 1983. Nonlocal stress field at Griffith crack, Cryst. Latt. Def. and 
Amorph. Mat. 10, 33-38. 

[8] Eringen, A.C., 1983. On differential equations of nonlocal elasticity and solutions of screw 
dislocation and surface waves. J. Appl. Phys. 54, 4703-4710. 

[9] Eringen, A.C., 1983. Interaction of a dislocation with a crack. J. Appl. Phys. 54, 6811-6817. 

[10] Eringen, E.C., 2002. Nonlocal Continuum Field Theories, Springer, New York. 

[11] Povstenko, Yu.Z., 1995. Circular dislocation loops in nonlocal elasticity. J. Phys. D: Appl. 
Phys 28, 105 111. 

[12] Povstenko, Yu.Z., 1995. Straight disclinations in nonlocal elasticity. Int. J. Engng. Sci. 33, 
575-582. 

[13] Povstenko, Yu.Z., Matkovskii, OA., 2000. Circular disclination loops in nonlocal elasticity. 
Int. J. Solids Structures 37, 6419-6432. 

[14] Lazar, M., 2003. A nonsingular solution of the edge dislocation in the gauge theory of 
dislocations. J. Phys. A: Math. Gen. 36, 1415-1437. 

[15] Kroner, E., 1963. On the physical reality of torque stresses in continuum mechanics. Int. 
J. Engng. Sci. 1, 261-278. 



26 



[16] Mindlin, R.D., 1964. Micro-structure in linear elasticity. Arch. Rat. Mech. Anal. 16, 51-78. 

[17] Mindlin, R.D., 1965. Second gradient of strain and surface-tension in linear elasticity. Int. 
J. Solids Structures 1, 417-438. 

[18] Mindlin, R.D., Eshel, N.N., 1968. On first strain gradients theory in linear elasticity. Int. 
J. Solids Structures 4, 109-124. 

[19] Green, A.E., Rivlin, R.S., 1964. Simple force and stress multipoles, Arch. Rat. Mcch. 
Anal. 16, 325-353. 

[20] Green, A.E., Rivlin, R.S., 1964. Multipolar continuum mechanics, Arch. Rat. Mech. 
Anal. 17, 113-147. 

[21] Gutkin, M.Yu., Aifantis, E.C., 1999. Dislocations in gradient elasticity. Scripta Mater. 40, 
559-566. 

[22] Gutkin, M.Yu., 2000. Nanoscopics of dislocations and disclinations in gradient elasticity. 
Rev. Adv. Mater. Sci. 1, 27-60. 

[23] Lazar, M., Maugin, G.A., 2003. Nonsingular stress and strain fields of dislocations and 
disclinations in first strain gradient elasticity. Int. J. Engng. Sci., in press. 

[24] Gutkin, M.Yu., Aifantis, E.C., 1996. Screw dislocation in gradient elasticity. Scripta 
Mater. 35, 1353-1358. 

[25] Gutkin, M.Yu., Aifantis, E.C., 1997. Edge dislocation in gradient elasticity. Scripta 
Mater. 36, 129-135. 

[26] Gutkin, M.Yu., Aifantis, E.C., 1999. Dislocations and disclinations in gradient elasticity. 
Phys. Stat. Sol. (b) 214, 245-284. 

[27] Aifantis, E.C., 2003. Update on a class of gradient theories. Mechanics of Materials 35, 
259-280. 

[28] Eringcn, A.C., 1984. Plane waves in nonlocal micropolar elasticity. Int. J. Engng. Sci. 22, 
1113-1121. 

[29] Povstenko, Yu.Z., 1998. Imperfections in non-local elasticity. J. Phys. IV France 8, Pr8- 
309-316. 

[30] Povstenko, Yu.Z., Matkovskii, OA., 1997. Edge dislocation in nonlocal elasticity medium 
with couple stresses. Mat. Met. Fiz.-Mekh. Polya. 40, 98-102 (in Ukrainian). 

[31] Lazar, M., Maugin, G.A., 2004. Defects in gradient micropolar elasticity: I. Screw disloca- 
tion. J. Mech. Phys. Solids 52, 2263-2284. 

[32] Lazar, M., Maugin, G.A., 2004. Defects in gradient micropolar elasticity: II. Edge disloca- 
tion and wedge disclination. J. Mech. Phys. Solids 52, 2285-2307. 

[33] Zubov, L.M., 1997. Nonlinear thoery of dislocations and disclinations in elastic bodies. 
Springer. 

[34] Altan, B.C., Aifantis, E.C., 1997. On some aspects in the special theory of gradient elas- 
ticity. J. Mech. Behav. Mater. 8, 231-282. 

[35] Kroner, E., 1981. Continuum Theory of Defects, in: Physics of Defects (Les Houches, 
Session 35), R. Balian et al., eds., North-Holland, Amsterdam, p. 215. 



27 



[36] DeWit, R., 1973. Theory of disclinations IV. J. Res. Nat. Bur. Stand. (U.S.) 77A, 607-658. 

[37] Lazar, M., 2003. Dislocations in the field theory of elastoplasticity. Comput. Mater. Sci. 28, 
419-428. 

[38] Eringen, A.C., 1999. Microcontinuum Field Theories I: Foundations and Solids, Springer, 
New York. 

[39] Nowacki, W., 1974. On discrete dislocations in micropolar elasticity. Arch. Mech. 26, 3-11. 

[40] Nowacki, W., 1986. Theory of Asymmetric Elasticity. Pergamon Press, Oxford. 

[41] Kessel, S., 1970. Spannungsfelder einer Schraubenversetzung und einer Stufenvcrsctzung 
im Cosserat Kontinmim. Z. Ang. Math. Mech. 50, 547-553. 

[42] Minagawa, S., 1977. Elastic fields of dislocations and disclinations in an isotropic micro- 
polar continuum. Appl. Eng. Sci. Lett. 5, 85-94. 

[43] Povstenko, Yu.Z., 1994. Stress functions for continua with couple stresses. J. Elast. 36, 
99-116. 

[44] Minagawa, S., 1979. Stress and couple-stress fields produced by circular dislocations in an 
isotropic elastic microploar continuum. Z. Ang. Math. Mech. 59, 307-315. 

[45] Kncsl, Z., Semela, F., 1972. The influence of couple-stresses on the elastic properties of an 
edge dislocation. Int. J. Engng. Sci. 10, 83-91. 

[46] Mindlin, R.D., 1963. Influence of couple-stresses on stress concentrations. Exper. Mech. 3, 
1-7. 

[47] Schaefer, H., 1962. Versuch einer Elastizitatstheorie des zweidimensionalen Cosserat- 
Kontinuums. Misz. Angew. Math., Festschrift Tollmien, Berlin, 277-292. 

[48] Carlson, D.E., 1966. Stress functions for plane problems with couple stresses. Z. Ang. Math. 
Phys. 17, 789-792. 



28 



